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Abstract 

We try to constrain the noncommutativity length scale of the theoretical model given in [l.) using 
the observational data from ACBAR, CBI and five year WMAP. The noncommutativity parameter 
is not constrained by WMAP data, however ACBAR and CBI data restrict the lower bound of its 
energy scale to be around 10 TeV. We also derive an expression for the amount of non-causality 
coming from spacetime noncommutativity for the fields of primordial scalar perturbations that are 
space-like separated. The amount of causality violation for these field fluctuations are direction 
dependent. 



f Electronic address: eakofor@phy.syr.edu 



^Electronic address: bal@phy.syr.edu 



^Electronic address: ajoseph@p hy.syr.edu| 



5 Electronic address: lppekows@phy.syr.edu 



^Electronic address: bqureshi@stp.dias.ie 



1 



I. INTRODUCTION 



In 1992, the Cosmic Background Explorer (COBE) satellite detected anisotropics in the 
CMB radiation, which led to the conclusion that the early universe was not smooth: there 
were small density perturbations in the photon-baryon fluid before they decoupled from 
each other. Quantum corrections to the inflaton field generate perturbations in the metric 
and these perturbations could have been carried over to the photon-baryon fluid as density 
perturbations. We then observe them today in the distribution of large scale structure and 
anisotropies in the CMB radiation. 

Inflation Q, H, fl stretches a region of Planck size into cosmological scales. So, at 
the end of inflation, physics at the Planck scale can leave its signature on cosmological scales 
too. Physics at the Planck scale is better described by models of quantum gravity or string 
theory. There are indications from considerations of either quantum gravity or string theory 
that spacetime is noncommutative with a length scale of the order of Planck length. CMB 
radiation, which consists of photons from the last scattering surface of the early universe can 
carry the signature of spacetime noncommutativity. With these ideas in mind, in this paper, 



we look for a constraint on the noncommutativity length scale from the WMAP5 



ACBAR 



12| and CBI 
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17l | observational data. 



In a noncommutative spacetime, the commutator of quantum fields at space-like sepa- 
rations does not in general vanish, leading to violation of causality. This type of violation 
of causality in the context of the fields for the primordial scalar perturbations is also dis- 
cussed in this paper. It is shown that the expression for the amount of causality violation is 
direction-dependent. 

In 181 ]. it was shown that causality violation coming from noncommutative spacetimes 
leads to violation of Lorentz invariance in certain scattering amplitudes. Measurements of 
these violations would be another way to put limits on the amount of spacetime noncommu- 
tativity. 

This paper is a sequel to an earlier work [l|. The latter explains the theoretical basis of 
the formulae used in this paper. In Q 

another approach of noncommutative inflation is 
considered based on target space noncommutativity of fields [2! 
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II. LIKELIHOOD ANALYSIS FOR NONCOMMUTATIVE CMB 



The CMBEasy 21] program calculates CMB power spectra based on a set of parameters 



and a cosmological model. It works by calculating the trans 
for scalar perturbations at the present conformal time rjo as 



er functions A; for multipole / 



22|] 



rvo 

Aj(fc, V = Vo) = I drj S(k, v)ji[ k (Vo ~ v)], 



o 



where S is a known "source" term and ji is the spherical Bessel function. (Here "scalar 
perturbations" mean the scalar part of the primordial metric fluctuations. Primordial met- 
ric fluctuations can be decomposed into scalar, vector and second rank tensor fluctuations 



according to their transformation properties under spatial rotations [23j. They evolve in- 
dependently in a linear theory. Scalar perturbations are most important as they couple to 
matter inhomogeneities. Vector perturbations are not important as they decay away in an 
expanding background cosmology. Tensor perturbations are less important than scalar ones, 
they do not couple to matter inhomogeneities at linear order. In the following discussion we 
denote the amplitudes of scalar and tensor perturbations by A s and At respectively.) The 
lower limit of the time integral in eq. (pQ) is taken as a time well into the radiation dominance 
epoch. Eq. (0Q) shows that for each mode k, the source term should be integrated over time 

7/. 

The transfer functions for scalar perturbations are then integrated over k to obtain the 
power spectrum for multipole moment /, 

Cf ] = (4tt 2 ) J dk k 2 P$ (k)\Mk,V = Vo)\ 2 , (2) 

where P$ is the initial power spectrum of scalar perturbations (cf. Ref. taken to be 

P$ (k) = A s k~ 3+ ^ ria ~ 1 ' with a spectral index n s . 

The coordinate functions x M on the noncommutative Moyal plane obey the commutation 
relations 

\x^x v \ = i9^ u , 9^ = -9 Ufl = const. (3) 

We set 0° = (6> 01 , 6> 02 , 9 03 ) to be in the third direction. In that case, 9° = 9 9° where the unit 
vector 9° is (0,0,1). 



We now write down eq. (79) of [lj], 

2 r 00 

(aima?w> 9 = ~ dk E ^'(-l) m (2/" + l)A; 2 A i (A;)A,(A;)P <I)0 (A;)v(^) 
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(4) 



where ij is the modified spherical Bessel function and H is the Hubble parameter during 
inflation. In the limit when 6 = eq. (jl]) leads to the usual CVs {24]: 

C < = ^tEM™>" = ( 47r2 ) / ^oWIAzCA;^ = m )\ 2 . (5) 

m 

Our goal is to compare theory with the observational data from WMAP5, ACBAR and 
CBI. These data sets are only available for the diagonal terms I — V of eq. and for the 
average over m for each I, so we consider only this case. Taking the average over m of eq. 
dlj), for lm = I'm' the sum collapses to 

cf = 2/TiX)< a ^>« = / dkk 2 P* (k)\MKv = Vo)W(0kH), (6) 

cf = a- (7) 

The CMBEasy integrator was modified to include the additional Iq code and the Monte 
Carlo Markov-chain (MCMC) facility of the program was used to find best-fit values for 9 H 
along with the other parameters of the standard ACDM cosmology. 

In the first run the parameters were fit using a joint likelihood derived from the WMAP5, 
ACBAR and CBI data. The outcome of this analysis was inconclusive, as the resulting value 
was unphysically large. This result can be understood by examining the WMAP5 data alone 
and considering a x 2 goodness-of-fit test, using 

2 \ ^ / ~~ ^l,data 



* 2 =£ 



07 



I 

where Ci^ata is the power spectrum and 07 is the standard deviation for each I as reported 
by WMAP observation. 
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FIG. 1: Transfer function A; for / = 10 as a function of k. It peaks around k = 0.001 Mpc . 

We expect noncommutativity to have a negligible effect on most of the parameters of the 
standard ACDM cosmology. We therefore consider the effect on the CMB power spectrum 
of varying only the new parameter HQ. To determine its effect, we consider the shape of the 
transfer functions A/ (k) as calculated by CMBEasy. The graphs of two such functions are 
shown in Figs. [T]and[2j As can be seen, these functions drop off rapidly with k, but extend 
to higher k with increasing /. (For example, in Fig. [U the transfer function for I = 10, 
A 10 , peaks around k = 0.001 Mpc -1 while in Fig. [2j the transfer function for I = 800, A 80 o, 
peaks around k = 0.06 Mpc -1 .) As io is a monotonically increasing function of k starting 
at io(0) = 1, this means that transfer functions of higher multipoles will feel the effect of 
noncommutivity first. 

The spectrum from the WMAP observation is shown in Fig. [3j Note in particular that 
the last data point, corresponding to I = 839 falls significantly above the theoretical curve. 
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FIG. 2: Transfer function A/ for I = 800 as a function of k. It peaks around k = 0.06 Mpc -1 . 

This means that x 2 can be lowered by a significant amount by using an unphysical value 
of HQ to fit this last point, so long as doing so does not also raise adjacent points too far 
outside their error bars. Performing the calculation shows that is indeed what happens. We 
therefore conclude that the WMAP data do not constrain HQ. 

Fig. H] shows the values of k which maximize Ai(k), as a function of I, which in turn gives 
a rough estimate of the region over which the transfer functions contribute the most to the 
integral in eq. (4), and hence the region over which changes in i (H8k) will most change the 
corresponding C\. Thus to improve the bound on HQ, we need data at higher I (I > 839). 
In addition, tighter error bars at these higher I will, of course, also help constrain the new 
parameter. 

Based on this analysis we performed a second run of CMBEasy excluding the WMAP 
data. This run resulted in a smaller, but still unphysically large, value of HQ. To see why this 
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FIG. 3: CMB power spectrum of ACDM model (solid curve) compared to the WMAP data (points 
with error bars). 



happens, we again consider the effect of varying only the new parameter HQ and examine 
the behavior of x 2 ■ 

ACBAR and CBI are CMB data on small-scales (ACBAR and CBI give CMB power 
spectrum for multipoles up to / = 2985 and I = 3500 respectively) and hence may be better 
suited to determination of HQ. A plot of x 2 versus HQ for ACBAR+CBI data is shown in 
Fig. [5j The plateau between HQ = Mpc and HQ = 0.01 Mpc is not physical, it results 
from limited numerical precision. Therefore, likelihoods calculated in this range only restrict 
HQ < 0.01 Mpc and hence cannot indicate whether the best fit is at HQ = Mpc or some 
small non-zero value. 

However, it is possible to put a constrant on the energy scale of spacetime noncommuta- 
tivity from HQ < 0.01 Mpc. We discuss this below. 
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FIG. 4: The values of k which maximize Aj(fc), as a function of I 

We can use the ACBAR+WMAP3 constraint on the amplitude of scalar power spectrum 
A s ~ 2.15 x 10~ 9 and the slow-roll parameter e < 0.043 [10] to find the Hubble parameter 
during inflation. The expression for the amplitude of the scalar power spectrum 



Tie VM 



where M p is the Planck mass, gives an upper limit on Hubble parameter: 



(9) 



H < 1.704 x 10 _5 M P . 



(10) 



On using this upper limit for H in the relation H9 < 0.01 Mpc, we have 9 < 1.84x 10 _9 m 2 . 

We are interested to know the noncommutativity parameter at the end of inflation. That 
is, we should know the value of the cosmological scale factor a when inflation ended. Most of 
the single field slow-roll inflation models work at an energy scale of 10 15 GeV or larger {24]. 
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FIG. 5: x 2 versus H9 for ACBAR data 



Assuming that the reheating temperature of the universe was close to the GUT e nerg y scale 



(10 16 GeV), we have for the scale factor at the end of inflation the value a ~ 10 -29 24|. Thus 
we have for the noncommutativity parameter, \J~Q < (1.84 a x lO -9 ) 1 / 2 = 1.36 x 10~ 19 m. 
This corresponds to a lower bound for the energy scale of 10 TeV. 



III. NON-CAUSALITY FROM NONCOMMUTATIVE FLUCTUATIONS 

In the noncommutative frame work, the expression for the two-point correlation function 
for the field tpe for the scalar metric perturbations contains hermitian and anti-hermitian 
parts [lj. Taking the hermitian part, we obtained the modified power spectrum 

P$ e (k) = P 9o (k) cosh( J ff^ t) • k), (11) 
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where P$ (k) is the power spectrum for the scalar metric perturbations in the commutative 
case (as discussed in H is the Hubble parameter during inflation. The constant spatial 
vector 9° is a measure of noncommutativity. The parameter 9 is related to 9° by 9° = 9z if 
we choose the z-axis in the direction of 6 , z being a unit vector. Also, 

-ikx 



$ fl (M) = J d 6 x<p e (x,t) e- lk ' x (12) 

This modified power spectrum was used to calculate the CMB angular power spectrum for 
the two-point tem per ature correlations. 

In this section 26], we discuss the imaginary part of the two-point correlation function 
for the field <pg. In position space, the imaginary part of the two-point correlation function 
is obtained from the "anti-symmetrization" (taking the anti-hermitian part) of the product 



of fields for a space-like separation: 

7;Mx,r/),<My,7/)]_ = ^(y^(x,f?)<My,?7) - w(y> v)M*, vi) • (13) 



The commutator of deformed fields, in general, is nonvanishing for space-like separations. 
This type of non-causality is an inherent property of noncommutative field theories con- 



structed on the Groenewold-Moyal spacetime 18] 



To study this non-causality, we consider two smeared fields localized at xi and X2. (The 
expression for non-causality diverges for conventional choices for P$ if we do not smear the 
fields. See after eq. (1201) .) We write down smeared fields at Xi and x 2 . 

^(a,x a ) = (^) 3/2 | d 3 x (p 9 ( X ) e - a ^\ (14) 

¥>(a,x 2 ) = (^) 3/ Y d 3 x ^(x) e~ a ^ 2 , (15) 

where a determines the amount of smearing of the fields. We have 

lim (®Y /2 f d 3 x ^(x) e -(x-x l} 2 = ^ (xi) . (16) 
a— >oo \7T/ J 

The scale 1 / yfa. can be thought of as the width of a wave packet which is a measure of the 
size of the spacetime region over which an experiment is performed. 
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We can now write down the uncertainty relation for the fields ip(oe, x x ) and ip(a, x 2 ) coming 
from eq. (|T3|) : 

A^(a, Xl )A^(a,x 2 ) > i| (0| [<p(a, Xl ), <p(a, x 2 )] |0) | (17) 

77ms equation is an expression for the violation of causality due to noncommutativity. 
We can connect the power spectrum for the field $o & t horizon crossing with the commu- 
tator of the fields given in eq. ( TT3l) : 



i<0[[^(k, 77),^(k',77)]_|0> 



(2vr) 3 P $0 (£;) smh(H(P ■ k) <5 3 (k + k'). 



horizon crossing 

Here we followed the same derivation given in [l|, using a commutator for the fields to start 
with, instead of an anticommutator of the fields, to obtain the above result. 
The right hand side of eq. (JTTj) can be calculated as follows: 



(0 1 [<f(a, xi),<p(a,x 2 



a\ 3 

71 J 

2 
2 



«\ 3 

7T 



cfarcfy <0|Mx),<My)]|0) e' 



-a(x-xi) 2 -a(y-x 2 ) 2 



d 3 £; <i 3 g 



(2tt) 3 

This gives for eq. (IT7j) 



^^I/T^tHs (0|[$ e (k),$ e (q)]|0) e-k-x-^-aKx-xO^Cy-x,) 2 ] 

^ J d 3 xd 3 yd 3 k P$ (k) smh(H(P • k) e -*-(x-y) e -a[(x- xl ) 2 +(y-x 2 ) 2 ] 
y rf 3 A; P$ (k) smh(H6° ■ k) e -^-*-(x 1 -x 2 )_ ^ 



A^(a,xi)Av?(a,x 2 ) > 



(27T) 



rf 3 A; P$ (A;) sinh(P^ • k) 



-ik-(xi— x 2 ) 



(20) 



The right hand side of eq. (1201) is divergent for conventional asymptotic behaviours of P$ 
(such as P$ vanishing for large k no faster than some inverse power of k) when a —>■ oo and 
thus the Gaussian width becomes zero. This is the reason for introducing smeared fields. 

Notice that the amount of causality violation given in eq. (1201 is direction-dependent. 

The uncertainty relation given in eq. fl20]) is purely due to spacetime noncommutativity 
as it vanishes for the case 9^ = 0. It is an expression of causality violation. 

This amount of causality violation may be expressed in terms of the CMB temperature 
fluctuation AT/T. We have the relation connecting the temperature fluctuation we observe 
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FIG. 6: The amount of causality violation with respect to the relative orientation between the 
vectors 6° and r = xi — X2. It is maximum when the angle between the two vectors is zero. Notice 
that the minima do not occur when the two vectors are orthogonal to each other. This plot is 
generated using the Cuba integrator [25f] . 



today and the primordial scalar perturbation $ 
AT(n,?7 ) ^ 

= = l^^mKWl 



rp / j -miy m , i-lmXP') 

Im 

— Mk^o^eimiik), (21) 

where n is the direction of incoming photons and the transfer functions A; take the primordial 
field perturbations to the present time r] . We can rewrite the commutator of the fields in 
terms of temperature fluctuations AT/T using eq. (j2ip . but the corresponding correlator 
differs from the one for the CMB temperature anisotropy. It is not encoded in the two-point 
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temperature correlation functions which as we have seen are given by the correlators of the 
anti-commutator of the fields. 

In Fig. [6l we show the dependence of the amount of non-causality on the relative orien- 
tation of the vectors Q° and r = x x — x 2 . The amount of causality violation is maximum 
when the two vectors are aligned. 



IV. CONCLUSIONS 



The power spectrum becomes direction dependent in the presence of spacetime non- 
commutativity, indicating a preferred direction in the universe. We tried a best-fit of the 
theoretical model in [l| with the WMAP data and saw that to improve the bound on HQ, 
we need data at higher I. (The last data point for WMAP is at I = 839.) We therefore 
conclude that the WMAP data do not constrain HO. We also see that tighter error bars 
at these higher / will also help constrain the noncommutativity parameter. The small-scale 
CMB data like ACBAR and CBI give the CMB power spectrum for larger multipoles and 
hence may be better suited for the determination of HQ. ACBAR+CBI data only restrict 
HQ to HQ < 0.01 Mpc and do not indicate whether the best fit is at HQ = Mpc or some 
small non-zero value. However, this restriction corresponds to a lower bound for the energy 
of Q of around 10 TeV. 

Further work is needed before rejecting the initial hypothesis that the other parameters 
of the ACDM cosmology are unaffected by noncommutivity. It requires performing a full 
MCMC study of all seven parameters. 

Also, we have shown the the existence and direction-dependence of non-causality coming 
from spacetime noncommutativity for the fields describing the primordial scalar perturba- 
tions when they are space-like separated. We see that the amount of causality violation is 
maximum when the two vectors, Q° and r = xx — x 2 , are aligned. Here r is the relative 
spatial coordinate of the fields at spatial locations xi and x 2 . 
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